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I.  INTRODUCTION 

This  report  describes  the  progress  made  during  the  period  of 
the  report  on  the  feasibility  study  and  design  of  an  Experimental 
Pulsed  Laser,  Remote  Cros.  wind  Measurement  System.  Work  during 
the  report  period  was  concentrated  on  the  analysis,  design  and 
construction  of  a sixty-four  element  optical  receiver,  described 
in  Sections  II  and  III,  and  an  analysis  of  the  accuracy  of  remote 
sensing  measurements,  described  in  Section  IV. 
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II.  RECEIVER  ANALYSIS  AND  DESIGN 
During  the  period  of  the  report  it  was  decided  to  replace  the 
silicon  vidicon  tube,  rotating  disk  receiver  described  in  our  last 
report  with  a photodiode  lens  array  receiver.  This  decision  was 
based  on  problems  encountered  with  the  silicon  vid icons. 

The  new  receiver  concept  has  many  advantages  including  the 
ability  to  work  off  of  a perfectly  diffuse  target,  no  moving  parts, 
much  closer  resemblance  to  an  actual  field  unit,  direct 
adaptability  to  a 1.06  urn  laser  wavelength  and  a much  lower  trans- 
mitted energy  requirement  for  a given  signal  to  noise  ratio. 

The  system  design  for  the  receiver  is  illustrated  in  Figures 
2-1  and  2-2.  The  heart  of  the  system  is  3 unit  cell  that  collects 
the  ruby  laser  energy  reflected  from  the  target  and  provides  an 
electrical  signal  suitable  for  storage  and  later  processing. 

Sixty-four  unit  cells  consisting  of  a plano-convex  lens,  two  baffle 
plates,  a photo  detector,  a current  to  voltage  transducer  and  an 
integrate  and  hold  amplifier  are  utilized  in  the  system.  The  lenses 
are  22.4  mm  in  diameter  with  a 150  mm  focal  length  and  are  arranged 
in  a two-dimensional  array  with  a center  to  center  spacing  of  30.48  mm. 

There  are  several  constraints  on  the  lens  spacing.  Since  we  are 
basically  measuring  the  covariance  function  with  the  instrument, 
the  spacing  should  be  such  that  three  or  four  spacings  span  the 
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majority  of  the  covariance  curve.  Also,  there  are  physical  con- 
straints on  the  total  array  size.  From  Figure  8 of  our  last  report, 
it  appears  that  a total  of  around  three  Fresnel  zones  are  required. 

The  lens  spacing  in  Fresnel  zones  of  course  depends  on  the  range  to 
the  target.  The  spacing  used  in  the  system  corresponds  to  0.917 
Fresnel  zones  at  1.6  Km  and  0.67  Fresnel  zones  at  3.0  Km. 

It  is  necessary  to  restrict  the  field  of  view  (FOV)  of  each 

unit  cell  in  order  to  limit  the  background  light  to  an  acceptable 

level.  Care  must  be  taken  however  to  not  overlv  restrict  the  FOV 

2 

or  fading  due  to  angle  of  arrival  fluctuation  may  introduce  errors 
into  the  measurements.  Fortunately,  neither  an  optical  filter  nor 
a precision  FOV  aperture  is  required.  The  lens  in  combination  with 
the  limited  size  of  the  active  area  on  the  detector  limits  the  FOV 
and  consequently  the  background  to  an  acceptable  level.  The  resultant 
FOV  is  approximately  six  milliradians. 

The  received  background  light  has  two  deleterious  effects.  It 
acts  as  a bias  that  could  saturate  the  unit  cell  if  not  restricted, 
and  its  fluctuations  constitute  a source  of  noise.  Each  of  these 
effects  is  evaluated  and  a signal  to  noise  ratio  analysis  for  the 
unit  cell  Is  contained  in  Section  III. 

The  photo  current  generated  in  the  detector  by  the  received 
laser  pulse  is  amplified  and  converted  to  a voltage  pulse  by  an 
operational  «mplifier  operating  in  the  current  mode  (refer  to  Figure 
2-2).  It  is  then  integrated  by  the  integrate  and  hold  amplifier  to 
provide  an  output  proportional  to  the  received  energy  in  the  pulse. 

Electronically  the  system  is  divided  into  four  independent 
modules,  each  containing  a four  by  four  array  of  unit  cells,  an 
analog  multiplexer,  a twelve  bit  analog  to  digital  converter  and  a 
temporary  storage  means  (refer  to  Figure  2-1).  The  outputs  of  each 
of  the  unit  cells  are  sequentially  multiplexed  during  the  hold  cycle 
to  the  analog  to  digital  converter  where  the  information  is  digitized 
and  placed  in  temporary  storage. 
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During  normal  operation  of  the  system,  a measurement  cycle 
starts  with  the  operator  relinquishing  control  to  the  microprocessor. 
The  microprocessor  then  charges  the  laser  power  supply  and  at  the 
aopropriate  times  provides  signals  to  fire  the  laser  flash  lamp  and 
Q switches,  to  integrate  and  hold  and  to  transfer  data.  Each 
nodule  goes  through  three  complete  cycles  for  each  measurement. 

The  first  cycle  occurs  before  the  laser  is  fired  and  collects  a 
set  of  data  (64  numbers)  representing  background  and  offset  in  the 
integrators.  The  second  and  third  cycles  occur  just  after  each  of 
the  two  laser  pulses  respectively.  After  the  end  of  the  third 
Ejdule  cycle,  the  data  in  temporary  storage  is  transferred  to  the 
microprocessor.  At  the  option  of  the  operator  it  can  be  punched 
on  paper  tape.  The  operator  also  has  the  option  of  whether  or  not 
to  subtract  the  background  and  offset  data.  The  option  was  provided 
so  that  during  alignment  and  test  the  subtraction  process  could  be 
bypassed. 

Timing  of  the  Q switch  pulses  and  the  integrate  cycles  is 
crystal  controlled.  The  delay  time  between  flash  lamp  excitation 
and  the  first  Q switch  pulse  is  adjustable  from  400  to  599  micro- 
seconds. The  timing  for  the  second  Q switch  pulse  is  designed  to 
provide  a fixed  delay  of  one  millisecond  after  the  first  Q switch 
pulse. 

An  analog  output  is  provided  to  display  the  data  collected  on 
an  oscilloscope.  It  can  be  used  to  preview  a data  set  before  deciding 
whether  or  not  it  should  be  punched  and  to  verify  that  signal  levels 
and  alignment  are  proper. 

An  input  connection  is  provided  to  allow  an  external  unit  cell 
to  be  used  for  recording  the  laser  output  energy  on  each  pulse. 

This  information  is  punched  directly  on  the  data  tape  and  greatly 
simplifies  calibration  of  the  system. 

The  finished  receiver  is  shown  in  Figures  2-3  through  2-6. 

The  electronics  fabrication  was  accomplished  utilizing  printed 
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Figure  2-4  Receiver , Front  View 
5 


circuit  boards  for  the  subassembly. 

In  addition  a master  printed  circuit  board  utilizing  an  array 
of  connectors  was  used  to  interconnect  the  various  subassembly 
boards.  This  greatly  reduces  the  amount  of  hand  wiring  and  in- 
creases the  reliability  ami  maintainability  of  the  system. 

Regulated  DC  power  and  other  interconnections  enter  the  receiver 
through  connectors  on  the  back  panel.  Additional  regulators  on  the 
circuit  boards  as  well  as  RFI  filters  on  the  back  panel  are  utilized 
to  minimize  interaction  between  the  receiver  and  the  laser  power 
supply. 

The  receiver  electronics  are  housed  in  an  aluminum  enclosure 
that  is  mounted  in  a cradle  with  two  degrees  of  freedom.  A telescope 
with  crosshairs  is  mounted  on  the  receiver  and  bore  sighted  with  the 
optics  to  allow  the  receiver  to  be  aimed  at  the  target.  A detailed 
description  of  the  receiver  circuitry  will  be  published  at  a later 
date. 
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III.  UNIT  CELL  ANALYSIS 

In  our  last  report,^"  it  was  noted  that  if  a gated  receiver  is 
utilized,  then  for  a given  pu3se  energy,  as  the  gate  time  is  de- 
creased, the  signal  to  noise  ratio  (S/N)  increases  until  signal 
shot  noise  limited  operation  is  achieved.  In  addition,  a crude 
estimate  of  the  S/N  versus  gate  time  was  developed.  A much  more 
exact  analysis  of  the  S/N  versus  gate  time  has  been  developed  for 
the  unit  cell  utilized  in  the  new  receiver  and  is  presented  in 
this  section. 

Figures  3-1  and  3-2  illustrate  the  electronics  approach  used 
in  the  unit  cell  design.  A photodiode,  operational  amplifier 
combination  operated  in  the  current  mode  is  utilized  as  a current 
to  voltage  transducer  (CVT) . The  voltage  output  of  this  combina- 
tion is  then  fed  through  a switch  to  an  operational  amplifier  con- 
nected as  an  integrator.  The  capacitor  used  for  integrating  has  a 
switch  across  it  that  is  used  to  "dump"  the  accumulated  charge 
after  a ir.  asureii._nt  has  been  made.  The  integrator  together  with 
its  switches  will  be  treated  as  a single  entity  that  acts  like  a 
"finite  time  integrator"  (FTI). 

The  random  process  problem  that  we  must  solve  is  illustrated 
in  Figure  3-3.  The  input  random  process  n(t)  with  autocorrelation 

function  R (t)  is  wide  sense  stationary,  "white"  and  represents 
nn 

any  of  the  noise  sources  associated  with  the  diode  or  the  CVT. 

The  impulse  response  hj (t)  represents  the  linea-  circuitry  between 

any  of  the  noise  sources  and  the  output  of  the  CVT.  The  exact 

forms  will  be  considered  later.  The  output  x(t)  is  then  wide  sense 

stationary  with  autocorrelation  function  R (t).  The  impulse  re- 

xx  2 

sponse  h2(t)  of  the  finite  time  integrator  is  given  by 

h2(t)  - G[u(t)  - u(t-T) ] (3-1) 

where  G is  a gain  constant  (1/RC)  and  u(t)  is  the  unit  step 
function.  In  general,  when  a stationary  process  is  sampled  or 
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switched,  the  output  becomes  non-st at ionary . Consequently, 

analysis  of  the  finite  time  integrator  must  be  approached  from 

the  standpoint  of  a non-stat ionary  process. 

For  a non-stat ionary  process,  the  autocorrelation  function 

4 

at  th2  output  of  a linear  system  is  given  by 


R (to ,t ! ) = R (tp , 1 1 ) * ho(t,) 

yy  ..  1 yX  '■  ‘ 


(3-2) 


where 

Ryx(t2,t1)  = Rxx(t2,t,)  * h2 (t i ) 

Working  first  with  (3-3)  and  utilizing  (3-1)  it  becomes 

a> 

RyX ( 1 2 , t ! ) = G f Rxx(t2-«)[u(ti-0t)  - u(t j -a-T) ]du 


where  use  has  been  made  of 
T = t2  ~ tj 

Making  the  change  of  variable 

P = t2  - a 

(3-4)  becomes 


R (t2 

yx  *■ 


t2-tl+T 


R (P)dP  = R (t2-t,) 

XX  yx  ^ 1 


(3-3) 


(3-4) 


(3-5) 
y* 

t2-tj 

From  (3-5)  it  is  concluded  that  the  cross  correlation  function 


RyX(t2,tj)  is  wide  sense  stationary. 

Starting  with  (3-2)  and  utilizing  the  results  of  (3-5)  and  the 
same  procedures  it  is  found  that 


R (t2,t,) 
yv  *■  1 


= G 


oo 

I 


(a-tj)[u(t2-a)  - u(t2-a-T)]da 


(3-6) 
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Making  the  change  of  variable 
Q = a - tj 


t2-tj-T 


Ryy(t2.tl>  = G 


J »yJtM><lQ  - By/t2-c,) 


t2-tl 


(3-7) 


Consequently,  the  output  of  the  finite  time  integrator  is  wide  sense 
stationary  which  greatly  simplifies  the  analysis. 

Since  the  output  y(t)  is  wide  sense  stationary,  the  variance 
can  be  expressed  as^ 


•if 


S (to)  did 

yy 


(3-8) 


where 


Syy(to)  = Snn(a»)!H1(to)|2|H2(to)|2 


(3-9) 


and  Hj(w)  and  H2(w)  are  the  Fourier  transforms  of  the  impulse  re- 
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sponses  h](t)  and  h2(t)  respectively.  For  the  finite  time  integrator 


|H2  (w)i 2 = G2 

(t)' 


(3-10) 


However,  |Hj(to)]2  depends  on  which  noise  source  is  being  evaluated. 
Assuming  that  each  noise  source  is  uncorrelated  they  may  be  analyzed 
separately  and  the  results  summed  to  give  the  total  variance  due  to 
the  noise  sources.  This  yields^ 


n 2t!R 


b f Hn2l*(")'2  ♦*gl»Ml»  + d“ 


(3-11) 


where 


R,‘ <*'  “ 
f n 


= (, .>?+..)  ?)? 


~ l+^Cf'Rp" 


;c(->  j-  = 


s ^[Rf/Rj+U' i fRi/(Rf+Ri)/((‘.-K'f)-  + H 


(">'%•] 2 


„ AlvT 

!  I  2 = — + 2e (I  +1  ,+It+ID) 
n rj  s d ■ B 


= Amplifier  voltage  noise  /V Hz 


r,  = Detector  shunt  resistance 
d 


R j *»  CVT  amplifier  input  resistance 


Rt  - Ri  rd/(Ri+rd> 


C.  = Total  capacitance  in  shunt  with  CVT  amplifier  Input 


»„  - [y<vci4C£»] 


if? 


. jj  = 2r.  . Gain-Bandwidth  product  of  the  amplifier 


I = Signal  Current 


I.  = Detector  dark  current 
d 


I = Detector  background  current 
15 


and 


I]  = CVT  amplifier  input  current 

The  parameters  R,  C and  R^.  are  defined  in  Figure  3-1. 


(3-12) 

(3-13) 

(3-14) 

(3-15) 

(3-16) 

terminals 

(3-17) 
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The  expression  given  in  (3-11)  neglects  any  noise  associated  with 
the  bias  source  and  assumes  that  C^  is  chosen  for  critically  damped 
operation. 

The  signal  to  noise  ratio  at  the  output  of  the  unit  cell  then 
is  given  by 


o 

0 * 

n 


(3-18) 


where  V is  the  response  of  the  unit  cell  to  the  signal.  Since 
out 

the  signal  current  pulse  is  very  short  in  comparison  to  the  rise 
time  of  the  amplifier,  it  behaves  like  an  impulse.  Consequently 
the  output  of  the  amplifier  is  just  the  impulse  response  of  the  net- 
work represented  by  (3-12).  The  output  is  given  by 


V(t)  = -URpRr1n?te“*nt 


u ( t ) 


(3-19) 


where  U is  the  total  received  energy  in  the  pulse  and  p is  the 
K 

diode  responsivity  in  amps  per  watt.  The  voltage  V(t)  is  integrated 
by  the  finite  time  integrator  to  yield 


out 


URpRf  ‘n 
RC 


L 


“,u  t 

t e n dt 


(3-20) 


For  T long  compared  to  1/m  , (3-20)  becomes 


out 


PURRf 

RC 


(3-21) 


In  order  to  calculate  the  signal  to  noise  ratio,  the  signal 
and  background  levels  must  be  evaluated.  The  received  signal  energy 
is  given  by 


U, 


U_T_ 
T R 


(3-22) 
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where 


U.j,  = Transmitted  energy 

T = Target  reflectance 

D^  = Receiver  diameter  in  meters 

and 

L = Range  to  target  in  meters 

Assuming  a power  incident  on  the  target  from  the  sun  of  0.03 
Watts/Cnr/urn,  the  background  power  level  is  given  by 

3 x 10fT  - 

PB Jf IF0V-0J2  mo  (3-23) 

where 

FOV  = D^/ Focal  Length 

BWq  = Optical  Bandwidth  in  ym 

and 

* Diameter  of  detector  active  area 

For  the  unit  cell,  the  optical  bandwidth  is  determined  by  the  spec- 
tral response  of  the  photodiode  and  is  approximately  0.65  urn.  Using 
a reflectivity  of  unity,  this  yields  for  the  unit  cell  a background 
power  level  of  approximately  0.7  microwatt.  For  a diode  responsivity 
of  0.4  amps/watt  and  a one  megohm  feedback  resistor  this  yields  .28 
u amp  for  I„  and  .28  volt  for  the  output  of  the  CVT. 

D 

The  DC  output  due  to  background  is  not  sufficient  to  saturate 
the  CVT,  but  would  cause  a relatively  large  offset  at  the  output  of 
the  integrator.  To  circumvent  this  problem,  AC  coupling  and  a "gate" 
switch  was  used  between  the  CVT  and  the  integrator.  In  picking  the 
low  frequency  cutoff  for  the  AC  coupling,  the  frequency  response  of 
the  CVT  and  signal  must  be  considered.  Figure  3-4  shows  the  impulse 
frequency  response  of  the  CVT.  As  can  be  seen  from  the  figure,  the 
low  frequencies  are  important  and  consequently  the  cutoff  for  the  AC 
coupling  should  be  made  as  low  as  is  practical.  For  the  unit  cell. 
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the  cutoff  is  approximately  2 KHz. 

The  power  spectral  density  at  the  output  of  the  CVT  and  the 
unit  cell  due  to  all  the  pertinent  noise  sources  has  been  calcu- 
lated in  accordance  with  (3-9)  through  (3-16)  and  is  shown  in 
Figure  3-5.  Utilizing  (3-18),  and  (3-20)  and  the  results  illus- 
trated in  Figure  3-5,  the  signal  to  noise  ratio  for  the  unit  cell 
has  been  calculated  and  is  shown  in  Figure  3-6  versus  gate  time 
for  a 0.1  reflectivity  target6  at  3 Km.  For  the  present  unit 
cell  design,  the  optimum  gate  time  is  around  six  microseconds 
which  yields  a signal  to  noise  ratio  of  approximately  53  db  at 
a range  of  3 Km.  Both  the  gate  delay  time  which  compensates  for 
signal  time  delay  due  to  target  range  and  the  actual  gate  time  are 
adjustable. 

In  the  present  unit  cell  design,  the  dominant  noise  sources 
are  the  CVT  amplifier  voltage  noise  and  background  cur- 
v ' rent  in  the  detector.  Further  improvement  in  the  signal  to  noise 

§,  *,  ratio  could  be  obtained  by  reducing  these  noise  sources  or  by 

utilizing  a shorter  gate  time.  The  latter  approach  would,  however, 
require  an  amplifier  with  a wider  bandwidth. 

The  unit  cell  design  was  tested  with  the  ruby  laser  and  good 
results  were  obtained.  An  oscilloscope  picture  of  the  output  of 
the  CVT  is  shown  in  Figure  3—7  for  a diffuse  target  at  500  meters. 
The  two  superimposed  traces  correspond  to  the  two  laser  output 
pulses. 


IV.  ACCURACY  OF  REMOTE  SENSING  MEASUREMENTS 

There  are  several  factors  associated  with  remote  sensing  util- 
izing a pulsed  laser  system  that  will  cause  errors  in  evaluating 
the  crosswind  velocity  and  need  to  be  investigated.  These  factors 
include  a partially  uncontrollable  source  si^e,  saturation  of  the 
scintillations,  receiver  calibration,  estimation  of  the  statistics 
of  the  received  intensity  utilizing  only  a finite  number  of  samples 
and  electronic  noise  in  the  system. 

As  discussed  in  our  last  report,  the  calibration  constant  that 
relates  the  measured  slope  of  the  covariance  function  to  the  path 
weighted  average  crosswind  velocity  is  a function  of  the  source  size 
at  the  target  end.  Unfortunately  the  source  size  is  a function  of 
the  turbulence  level  and  the  degree  of  coherence  of  the  laser  trans- 
mitter. The  coherence  of  the  laser  transmitter  will  vary  in  an  un- 
controllable manner  from  measurement  to  measurement  and  from  pulse 
to  pulse  in  a single  measurement.  Preliminary  analysis  of  the  effect 
on  the  calibration  constant  indicates  that  as  the  source  size  in- 
creases, the  sensitivity  of  the  calibration  constant  to  a changing 
source  size  decreases.  This  is  illustrated  in  Figures  4-1  through 
4-3.  By  proper  choice  of  transmitter  beam  divergence  and/or  by 
estimating  the  source  size  from  the  measured  data  it  should  be 
possible  to  reduce  the  effect  on  crosswind  measurement  accuracy  to 
a reasonable  level.  At  this  point  the  work  has  not  proceeded  far 
enough  to  warrant  presenting  a detailed  analysis.  This  will  be  done 
in  the  fim.1  report  for  the  year. 

Saturation  of  the  scintillations  will  cause  the  calibration 
constant  to  change  with  turbulence  level.  Work  by  NOAA  indicates 
that  a proper  choice  of  source  and  receiver  sizes  will  minimize  the 
effect.  In  addition,  by  estimating  the  degree  of  saturation  from 
the  measured  data,  it  may  be  possible  to  compensate  for  the  effects 
of  saturation.  Considerable  additional  work  is  required  on  this 
problem. 
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Each  cell  in  the  receiver  array  must  be  calibrated  to  compen- 
sate for  variation  of  sensitivity  from  cell  to  cell.  This  is 
accomplished  by  utilizing  a diffuse,  poorly  reflecting  target  at 
a range  of  about  500  feet.  This  provides  a very  nearly  uniform 
illumination  level  to  the  receiver  for  calibration  with  a minimum 
of  turbulence  effects.  Multiple  sets  of  data  are  taken  and  pro- 
cessed on  a digital  computer  to  obtain  calibration  factors  accurate 
to  within  a few  percent. 

Remote  sensing  by  utilizing  the  effect  of  atmospheric  turbu- 
lence on  a propagating  wave  requires  that  certain  statistical 
parameters  be  measured.  For  the  specific  case  of  crosswind 
measurement  these  parameters  are  the  mean,  variance,  covariance 
function  and  slope  of  the  covariance  function.  Due  to  constraints 
on  the  time  allowed  for  the  measurement  and/or  the  number  of 
spatial  locations  that  can  be  utilized,  these  parameters  must  be 
estimated  by  utilizing  a finite  number  of  spatial  and/or  temporal 
samples.  Consequently,  this  limits  the  accuracy  with  which  the 
crosswind  may  be  evaluated  and  in  fact  may  be  the  limiting  effect 
with  respect  to  random  errors. 

In  order  to  analyze  the  error  due  to  using  a finite  number  of 
samples,  higher  order  joint  moments  of  the  intensity  must  be.  evalu- 
ated. Assuming  that  the  joint  distribution  for  the  intensity  is 
log-normal,  the  higher  order  moments  may  be  evaluated  in  terms  of  the 
mean,  variance  and  covariance.  This  is  fortunate  in  that  the  high- 
er order  joint  moments  for  the  intensity  have  not  in  general  been 
evaluated  theoretically. 

The  joint  moments  to  be  evaluated  are  of  the  form 
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where  Pj(«»)  is  the  Nth  order  joint  density  function  for  the  intensity 
and 
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using 


(4-2) 


this  can  be  expressed  as 
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where  is  the  nultivariate  jointly  noraal  density  function 
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Combining  these,  a rather  general  nth  moment  formulation  for  the 
intensity  is  obtained. 
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The  evaluation  of  this  integral 
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which  may  be  written  in  a slightly  more  usable  form  as 
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This  formulation  for  higher  order  moments  of  the  intensity  is  ex- 
pressed in  terms  of  first  and  second  moments  of  the  log-intensity. 
With  suitable  manipulation  it  can  however  be  expressed  in  terms  of 
first  and  second  moments  of  the  intensity.  Utilizing  (4-10)  to 
evaluate  the  first  and  second  moments  of  the  intensity  yields 
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which  can  be  manipulated  to  generate 
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Utilizing  (4-14)  and  (4-15)  in  equation  (4-10)  yields 
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Equation  (4-16)  can  be  used  to  evaluate  an  Nth  order  joint  moment  of 
the  intensity  utilizing  only  the  means,  variances  and  covariances  of 
the  intensity.  For  the  case  where 

<ii2>  = <i  2> 


n n( 1 « 

n n l — + l 

i=*j  \<I>2 


(4-17) 


Formulations  for  the  sampled  mean,  variance,  covariance  function 
and  slope  of  the  covariance  function  contained  in  our  last  report*  are 
repeated  here  in  a slightly  modified  form  for  convenience. 
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An  analysis  of  the  root-mean-square  (RMS)  error  due  to  utilizing  a 
finite  number  of  samples  to  estimate  the  statistical  parameters 
associated  with  the  turbulence  has  been  performed  for  the  sample 
mean  and  sample  variance.  It  is  assumed  that  an  n x n two-dimensional 
spatial  or  space-time  array  is  used  to  collect  the  sampled  values. 

The  mean  square  error  for  the  sample  mean  can  be  expressed  as 
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Using  equation  (4-18)  for  the  sample  mean 
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which  can  be  manipulated  to  yield 
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Again  using  (4-22) 
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Combining  (4-24)  and  (4-25)  the  mean  square  error  for  the  sample 
mean  becomes 
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where  the  covariance  is  given  by 
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For  the  case  of  a spatial  array  with  equal  normalized  spacing  S* 
between  cells  this  becomes 
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(4-28) 


Under  this  condition  and  after  considerable  manipulation  the  mean 
square  error  normalized  to  <I>2  becomes 
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where  C is  the  normalized  covariance.  The  results  of  applying 
N 

2 

(4-29)  to  the  case  of  a spherical  wave  are  illustrated  in 
Figures  4-4  and  4-5.  As  can  be  seen  from  the  figures, 
the  error  due  to  using  a finite  number  of  samples  can  be  consider- 
able. From  (4-2?;  it  can  be  seen  that  as  the  variance  of  the 
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intensity  increases  and  as  the  covariance  curve  spreads  our.,  t he 
error  for  a given  array  size  increases.  Consequently,  as  the 
source  size  increases,  the  variance  decreases  lowering  the  error; 

i 

however,  the  covariance  curve  spreads  out,-  thereby  tending  to  in- 
crease the  error.  The  covariance  curve  also  spreads  out^*"  with  in- 
creasing levels  of  saturation. 

Equation  (4-29)  can  also  be  used  to  evaluate  the  effect  of 
electronic  noise  on  the  error.  The  electronic  noise  is  composed 
of  thermal  noise,  and  Poisson  fluctuations  of  bias  current  and 
signal  current.  The  thermal  and  bias  current  noise  is  independent 
of  the  signal  and  its  turbulence  induced  fluctuation  and  not 
correlated  from  cell  to  cell.  This,  however,  is  not  the  case  for 
the  shot  noise  associated  with  the  signal.  In  order  to  simplify 
the  analysis,  it  will  however  be  assumed  that  the  electronic 
noise  and  turbulence  fluctuations  are  not  correlated. 

With  this  assumption. 


and 
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where  <?e2  is  the  variance  due  to  electronic  noise.  This  yields  from 
(4-26) 


<I>2  <I>2n2  <I>2  n2  S/N  <I>2 


(4-30) 


For  this  case,  the  errors  due  to  the  turbulence  fluctuations  and 
electronic  noise  add  in  the  mean  square  sense.  As  an  example  and 
considering  only  the  electronic  noise,  a signal  to  noise  ratio  of 
only  20  db  and  a 64  element  array  yields  a normalized  error 
due  to  electronic  noise  of  only  1.25% 

The  normalized  mean  square  error  for  the  sample  variance  can 
be  expressed  as 
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Using  equation  4-19  for  the  sample  variance, 
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The  higher  order  moments  of  the  intensity  included  in  (4-36)  can 
be  evaluated  utilizing  (4-17).  They  are  given  by 
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Equation  (4-31)  in  conjunction  with  (4-36)  through  (4-40)  gives 

the  normalized  mean  square  error  for  the  sample  variance.  The 

2 

results  of  applying  this  to  the  case  of  a spherical  wave  is  illus- 
trated in  Figures  4-6  and  4-7.  It  is  noted  that  the  normalized 
sample  error  increases  as  the  variance  of  the  intensity  increases 
and  as  the  covariance  increases.  Therefore,  the  general  discussion 
on  the  effect  of  source  size,  turbulence  level  and  saturation  on 
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the  sample  mean  error  also  applies  to  the  sample  variance  error. 

The  effect  of  electronic  noise  on  the  sample  variance  error 
can  be  evaluated  using  (4-31),  (4-32),  and  (4-36).  Again  assuming 
no  correlation  between  the  turbulence  fluctuations  and  the  electronic 
noise 
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The  errors  due  to  the  turbulence  fluctuations  and  electronic  noise 
add  in  the  mean  square  sense.  The  additional,  normalized  mean  square 
error  due  to  electronic  noise  is  given  by 
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As  might  be  expected,  the  electronic  contribution  to  the  error  in- 
creases as  Oj2  decreases.  Conversely  the  error  due  to  turbulence 
fluctuations  increases  as  increases  ar.d  the  two  effects  compete. 
It  is  expected  however  that  at  reasonable  electronic  signal  to  noise 
ratios,  the  error  due  to  turbulence  fluctuations  will  dominate. 
Equation  (4-43)  can  be  simplified  and  put  in  terms  of  the  electronic 
signal  to  noise  ratio  as  follows: 


32 


A • I 1 - /.  l\  <l  :'  °T?  # J_ 

ir  -’j*  S/S  y n-y  -Jj2  j2  * S/N 


^ + £)  ir 


(S/N)2 


(4-44) 


The  contribution  of  electronic  noise  to  the  normalized  variance 
error  is  illustrated  in  Figures  4-8  and  4-9. 

The  analysis  of  errors  associated  with  remote  sensing  of  the 
crosswind  is  not  yet  completed.  Remaining  work  includes  formulating 
the  same  covariance  error,  sample  slope  of  the  covariance  error 
and  crosswind  measurement  error  and  applying  it  to  the  actual  pulsed 
laser  remote  crosswind  sensor  configuration. 
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